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ABSTRACT 

In this publication electroweak next-to-leading order corrections to semileptonic £>-meson decays 
into (pseudo) scalar final states are presented. To this end, these corrections of O(aGp) have 
been calculated in the QED-enhanced phenomenological model, incorporating the bound-state 
mesons as its degrees of freedom, and matched to a similar calculation on the level of constituent 
partons in the full Standard Model. Consequently, the effects arising due to corrections of the 
partial decay widths on the extraction of the CKM matrix elements |V^b| an d | Vub| are detailed. 
Further, the results of two independent Monte-Carlo implementations are presented: One is 
the dedicated, strict fixed-order generator BLOR, and the other is embedded into the generic 
Yennie-Frautschi-Suura-type resummation of PHOTONS++, which is part of the multi-purpose 
event generator SHERPA. The resulting distributions are compared against the standard tool 
used in many experimental analyses, PHOTOS, showing improvements on the shapes of kinematic 
distributions of both the lepton and the final state meson. 
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1. Introduction 

In the Standard Model, the Cabibbo-Kobayashi-Maskawa (CKM) matrix [UI2] governs the 
charged current weak interactions between the up- and down-type quarks of the three fermion 
generations. The precision determination of its matrix elements and its CP-violating complex 
phase in the B meson sector has been the focus of intense research over the past decade. The 
combination of various measurements to test the unitarity of the CKM matrix is considered a 
strong instrument in the search for physics beyond the Standard Model OH]. 
In the present paper, a calculation of the electroweak next-to- leading order corrections in ex- 
clusive semileptonic B meson decays into (pseudo) scalar mesons, B — )■ D £ B — > Dq£u£ and 
B — > tt£i'£, £ denoting either an electron or a muon, is presented. Next-to-leading order cor- 
rections to such decays are an important aspect in the extraction of the CKM matrix elements 
|V^b| an d |Vub| at .B-factory experiments. Virtual electroweak bosons running in the loop as 
well as real photon emissions off all charged particles present in the decay alter the resulting 
decay dynamics and enhance the weak decay rate. To correct for the changed decay dynamics, 
experimentalists use approximative all-purpose next-to-leading order algorithms. These exploit 
universal factorisation theorems in the soft and/or collinear photon energy limit [5HZ]- hi addi- 
tion, the total hadronic decay rate of semileptonic decays are corrected by the known leading 
logarithm of the virtual corrections of the partonic decay [HE] • 
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Experience from exclusive semileptonic K meson decays illustrate the importance of having a 
good understanding of such radiative effects: until 2004 the global average of the extracted 
value of |Vus| from and Kf^ decays implied the violation of CKM unitarity by two standard 
deviations [10]. Further measurements proved dissonant with these findings [TlHl4] . indicat- 
ing that the achieved experimental precision needed an improved understanding of electroweak 
corrections. Since for many decays next-to- leading order calculations do not exist, experiments 
often use the approximative all-purpose algorithm PHOTOS |15t[16] to study the reconstruction 
efficiency and acceptance. The accuracy of this approach was tested by the KTeV collabora- 
tion, using the measured photon spectra from radiative Kf ?j decays: the angular distribution 
of the simulated photons did not agree well with the predicted spectrum [TTJ. This lead to 
the development of the next-to-leading order Monte Carlo generator KLOR (see [E]), whose 
next-to-leading order calculation is based on a phenomenological model. Its predicted angular 
photon distribution agreed satisfactorily with the measured spectra. Although this approach 
is very precise, it is also the most complicated one to adopt for an experiment: electroweak 
next-to- leading order calculations only exist for a few decay modes, sometimes only valid in a 
limited region of phase-space. Most of these calculations are evaluated numerically and rely on 
customised Monte Carlo generators. 

Over the last 10 years, an increasing amount of data and a better understanding of detector ef- 
fects lead to a very accurate picture of physics at the £>-factory experiments. This increased pre- 
cision then lead to the demand of knowledge of next-to-leading order electroweak effects beyond 
the precision of approximative all-purpose algorithms. The present paper aims at improving the 
status quo by providing a prediction for both total decay rates and differential distributions of a 
few representative kinematic variables. For the latter the predictions of a dedicated Monte-Carlo 
generator, BLOR [18j, is compared to two all-purpose generators, SHERPA/PHOTONS++ pT9| [20] 
and PHOTOS |15lll6j. While the latter is a QED-parton-shower Monte Carlo program intended 
to supplement generic leading logarithmic corrections to pure leading order decay generators, 
the former is a full-fledged hadron-level Monte Carlo generator for collider physics whose inter- 
nal leading order (hadronic) decays are supplemented by a universal soft-photon-resummation 
systematically improved, where possible, by known exact next-to-leading order matrix elements. 
While the improved description of inclusive decay rates directly gives small corrections to the 
extracted values of \V C ^\ and |y u b| from semileptonic decays, the improved description of the 
decay kinematics influence extrapolation to corners of the phase space and, therefore, leads to 
both direct and indirect corrections. 

The considerations of the present paper proceed as follows: Sec. 12.11 briefly reviews exclusive 
B — > X£u£ decays at tree-level. Thereafter, Sec. 12.21 develops the next-to-leading order for- 
malism, reviewing both the partonic short-distance results of [E1I2], the hadronic long-distance 
QED-improved effective decay and their matching to one-another, including also a detailed dis- 
cussion on non-universal structure-dependent terms in Sec. 12.31 This model is then embedded 
into the resummation in the soft limit of Yennie, Frautschi and Suura [6] in Sec. 12.41 Sec. [3] then 
shortly reviews the basic principles of both BLOR and SHERPA/PHOTONS++ where the calcula- 
tions of Sec. [2] have been implemented, and of PHOTOS. The total inclusive decay rates obtained 
are shown in Sec. 14.11 while differential distributions are shown in Sec. 14.21 a ls° detailing the 
improvement over the current estimates. The influence of the structure-dependent terms, where 
known, on the results is presented in Sec. 14.31 Sec. [5] finally summarises the results. 
Note that the charge-conjugated modes are implied throughout the present paper. 
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Fig. 1 The tree-level weak B — > X £ + v% decay is shown both in the phenomenological picture (left) 
and, at parton level, in Fermi's theory as low energy approximation of the Standard Model 
(right). The shaded circle represents the effective vertex parametrised by form factors f±, 
x 6 {u, c}. 



2. Phenomenological model 

2.1. Tree-level decay revised 

The phenomenological interaction Lagrangian of the weak B — )■ X tv decay to a (pseudo) scalar 
final state in Fermi's theory, with constant form factors of the hadronic current, i±, is given by 

C w = ^|Kb [(f+ + l-)<i>x d^ B + (f+ - f-)0B d^x] $ v Panvrtt + h.c. , (2.1) 

where ipg and t/j u are the Dirac fields of the lepton and the neutrino, <f>s and (px are the scalar 
fields of the initial and final state mesons, Gp the Fermi coupling, the CKM matrix element 
governing the strength of the b — > x transition, and P& = I + 75 is derived from the right-handed 
projection operator by absorbing the factor i into the coupling definition. The Lagrangian of 
eq. (|2.ip leads to the transition matrix element 

M° = -i-jLv A Hp(pB,Px\t) u u P R ^v e , (2.2) 

with the hadronic current, generalised to variable form factors, 

H^ PB ,Px;t) = {X\^ x P Kl ^ b \B) = ( PB +px) fM U(t) + {PB-px)^-(t). (2.3) 

The four-momenta labels in eq. (j2.2[) are introduced in Fig. [TJ The generalised form factors 
f-t = f±(t) now describe the phase-space dependent influence of the strong interaction on the 
weak decay dynamics and are functions of the squared momentum transfer from the hadronic 
to the leptonic system only, given at tree-level by 

t = {PB-Pxf = {pi+Pvf ■ (2.4) 
The tree-level differential decay rate in the B-meson rest frame is then given by 

dr o = \M° \ 2 dE x dE e , (2.5) 

with Ex = Px and Ei =p\. The explicit expressions of the f±(t) as function of the momentum 
transfer squared for the processes considered in this paper, B — > D lv, B — > D^lu and B — > 
7r Iv, can be found in App. fA]. 



throughout this paper denote a matrix element at 0(Gf a n ) with m photons in the final state. The total 
decay rate at 0(Gf « n ) is denoted as T 7 ^. 
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2.2. Next-to-leading order corrections 



The arising electroweak next-to-leading order corrections can be divided into two energy regimes: 
short-distance corrections at parton level, and long-distance corrections within the phenomeno- 
logical model. First, Sec. 12.2.11 will discuss how both descriptions can be matched and renor- 
malised. Sec. 12.2.21 then reviews the calculation of the virtual short-distance corrections of [BUS]- 
The long-distance corrections, following from an extension of the phenomenological model, are 
then discussed in Sec. 12.2.31 



2.2.1. Matching of different energy regimes 

The aim of this section is to develop a formalism to calculate the corrections at 0(aG-p). The 
standard approach involves calculating the one-loop graphs for the B — > X I vi decay in the 
effective theory with counterterms and compare it to the renormalised Standard Model result. 
Fixing the counterterms results in the desired matching of both results. The effective theory 
itself is non-renormalisable, but the Standard Model can be renormalised to measured quantities, 
e.g. the Fermi coupling constant of the muon decay, the electron mass, and the fine-structure 
constant, in order to produce finite predictions. Such a matching procedure was carried out 
in great detail by [2T] for semileptonic Kaon decays where the leading order phenomenological 
decay is described by a chiral Lagrangian. 

In the present paper, however, an alternative route is pursued. Consider a general logarithmi- 
cally divergent iV-point tensor integral of rank p with a single massless photon propagator. It 
can be cast in the form 

/j-mi . . . k^p 
* k V3T^: - (2 ' 6) 

with denominators di = (pi — k) 2 — m|. The integral can then be split according to 

fcw • • • k^ 1 ■ ■ ■ kf 1 ? 



r W "*(pi PAT-l) OC J d 4 k 



k 2 di--- d7v_i [k 2 - A 2 ] di ■ ■ ■ d N . 



■ ■ ■ h^P 

d 4 fc r?2 * - . (2.7) 

[k 2 - A 2 } di • • • d N -x 

This amounts to regulating the ultraviolet behaviour of the first term using an unphysical 
photon-like vector field of mass A and opposite norm, as proposed by Pauli and Villars in |22| . 
Its infrared behaviour is left unchanged, thus, relying on the Kinoshita-Lee-Nauenberg theorem 
[231124] . these divergences are left to be canceled by the real corrections. The second term of 
eq. (|2.7p is the equivalent of eq. ()2.6p . this time only with a massive photon. Hence, it is infrared 
finite and possesses the identical ultraviolet behaviour. 

Transferring this observation to the present case of semileptonic B meson decays where, both in 
the effective theory and in the Standard Model, there is at most one massless photon propagator 
in any one-loop diagram, the virtual emission matrix element can be decomposed as 

Ml = Ml M (A)+Ml sd (A). (2.8) 

The term Ai^ ld is now comprised of the Pauli- Villars regulated exchange of a massless photon, 
including its infrared divergence. The specific UV regulator effectively restricts the virtual 
photon's momentum to be smaller than A. Hence, it describes long-distance (ld) interactions 
only. 

The term A^J sd , on the other hand, carries the full ultraviolet behaviour of Mq- It thus can be 
used for renormalising all parameters. Consequently, because eq. ()2.7p is exact, all parameters in 
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Mq ld are then renormalised automatically. Through the photon mass, its virtual propagator's 
momentum is effectively restricted to be larger than A. Hence, this term describes the short- 
distance (sd) interactions only. 

The above is exact as long as the same Lagrangian input is used to calculate both the short- 
distance and the long-distance parts. In practice, however, due to the confining, non-perturbative 
nature of Qcd this is not feasible for the processes at hand. For scales larger than the hadron 
mass, its parton content can be resolved, electroweak corrections have to be calculated on the 
basis of (constituent) quarks. For scales smaller than the hadron mass, its parton content 
cannot be resolved, the bound-state hadrons themselves are the relevant degrees of freedom. 
Thence, supposing A is set such, that it effectively separates those two regimes, the long-distance 
Qed corrections .Mq ld can be calculated using the phenomenological model and for the short- 
distance corrections .Mq sd the full Standard Model has to be invoked. This is justified, in 
principle, by the assumption, that the phenomenological model describes the Standard Model 
and its effective degrees of freedom at these low scales. This mere fact, however, directly leads 
to inconsistencies at the matching scale A, where both models should give the same answer. 
Thus, this matching is only approximate and the associated systematic uncertainties to this 
method can be estimated by varying the matching parameter, cf. Sec. 14.11 
The above reasoning leads to an optimal value for A: the smallest hadronic mass in the decay. 
Then, only low-energy virtual photons, not able to resolve either of both mesons, are described 
by the effective theory, while high-energy virtual photons are described by the short-distance 
picture of the full Standard Model, resolving the partonic content of both the charged and 
the neutral meson involved. Further, as long as A > -E™ ax , the kinematic limit of the photon 
energy in single photon emissiorH, the real emission of photons off these charged mesons are 
also correctly described by the phenomenological model (except for structure dependent terms 
discussed in Sec. 12 .3[) . 

Nonetheless, it has to be noted that there are conceptual problems if both hadronic scales 
differ significantly. Then, there is a large intermediate regime, where virtual photons are able 
to resolve one meson, but not the other. By the above choice of A it is expected to give the 
best approximate description in this region. Further, if the third scale £™ ax exceeds A, real 
radiation, in the present ansatz always described using the phenomenological model, is able to 
resolve the final state meson, be it charged or neutral, as well. However, even if a considerable 
fraction is radiated at scales above A, this should have only negligible effects on the total decay 
rate as the bulk of the radiation is in the region k — > and, therefore, adequately described. 

2.2.2. Short-distance next-to-leading order corrections 

The well-known Standard Model Lagrangian is used in this study to describe the partonic 
b — > xtv decay in the short-distance regime. A representative collection of relevant next-to- 
leading order corrections to the tree-level decay, involving the exchange of virtual photons, W 
and Z bosons as well as Higgs scalars, is depicted in Fig. [2j Besides vertex corrections to the 
b-x-W and l-v-W vertices, wave function and propagator corrections, box diagrams involving 
the additional exchange of a neutral 7, Z or h bosons between the hadronic and leptonic systems 
are present. These next-to- leading order corrections read, calculated in [8] within the current 
algebra framework and concentrating on the renormalisation of the bare Fermi coupling G-p for 



2 The maximum photon energy for single photon radiation in the rest frame of a decaying particle is half its 
mass, neglecting all other decay products' masses. Allowing for massive decay products further reduces this 
kinematic limit. 
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Fig. 2 Representative Feynman diagrams of the Standard Model partonic decay b — > x£ + vi are 
shown. The white circles indicate hadronic contributions that are neglected in the short- 
distance expansion. 



this process, to leading logarithmic accuracy 
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(2.9) 



with A4g = GfMq, i.e. the leading order matrix element stripped of the Fermi coupling constant, 
and Q being the average charge of the quark line. The ellipsis stands for non-logarithmic terms. 
For the photonic contributions an infrared regulator has been introduced in the form of a 
photon mass A <C mw- Thus, the resulting matrix element has exactly the form required for 
the matching outlined in Sec. 12.2.11 The arising loop corrections are ultraviolet (UV) divergent 
and have been regularised in [8] by a UV cutoff set to mw- Renormalisation of the parameters, 
again focusing on , is then achieved by comparison to the muon decay computed in the same 
computational framework, yielding the relation 
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where Gf is now the renormalised Fermi decay constant as measured in the muon decay. This 
leads to the partonic short-distance virtual matrix element in the renormalised theory 
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(2.11) 



where Gf has been reabsorbed into the leading order matrix element -M[j. In the case of 
semileptonic B decays Q = \\Qi + Q&\ = g, x G {u, c}, this gives 



- In — • Mn + .... 

7T A U 



(2.12) 



The logarithm in eq. (|2.12p then represents the leading logarithmic corrections to 0(ckGf) due 
to virtual particle exchange with (virtual) photon energies above A. 
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Fig. 3 The Feynman diagrams for the next-to-leading order corrections to B + — > X° ^ + v decays are 
shown. 



2.2.3. Long-distance next-to-leading order corrections 

The Qed long-distance corrections to the phenomenological hadron decay can be calculated in 
an effective model that arises by requiring the phenomenological Lagrangian of the leading order 
decay to be invariant under local U(l) em gauge transformations. Assigning the usual charges 
the following interaction terms in the Lagrangian arise in addition to eq. (|2.1p 

/Wd = - eQfWI>tAp - ieQ^A^+d^- - ^8^+) + e % Q\A^<^^ 

+ ieV2G F V xy t±(Q B ±Q x )<pB<PxA^ u P R ^ e + h.c. , (2.13) 

wherein the summation over (f) £ {<j>B,<l>x} is implied. In addition to the point-like lepton- 
photon and meson-photon interactions, a vertex emission term arises. This term is connected 
to the bound-state nature of the meson. It is infrared finite and needed for gauge invariance. 
Further, in eq. f|2. 13j) it is assumed that the meson-photon interaction is sufficiently described 
by scalar Qed. Additional terms arise when moving away from this assumption, including 
intermediate lines of excited hadrons necessitating X* — > vertices as well as contributions 
due to off-shell currents. These terms are discussed on general grounds in Sec. 12.31 and will 
be largely neglected in this study. This, in most cases, roots in their unavailability or, where 
known, in their smallness. 

Figs. [3] and H] depict the relevant real and virtual diagrams for B + — > X° £ + v and B° — > X~ i + v 
decays at 0(a Gp). The real corrections diagrams a, b and c correspond to the emission of a real 
photon from either the charged legs of the decay, or the charged vertex itself. The virtual correc- 
tions group into three categories: diagrams d and e concern the wave- function renormalisation 
of the charged legs while diagram f is the dominant inter-particle photon exchange. Diagrams 
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Fig. 4 The Feynman diagrams for the next-to-leading order corrections to B° —> X £ + v decays are 
shown. 



g and h are again due to emissions off the charged effective vertex and are, thus, infrared finite. 

The corresponding subtleties involving the vertex emissions are detailed in App. |Bj 

In the virtual amplitude A4q ld the arising ultraviolet divergences are regularised using the Pauli- 

Villars prescription [22] by introducing an unphysical heavy photon of mass A and opposite 

norm. Consequently, the virtual corrections have exactly the form required by the matching 

procedure outlined in Sec. 12.2. II and all real emission processes are described in the long-distance 

picture. 

Squaring the phenomenological real and virtual matrix elements results in the real and virtual 
next-to-leading order differential rates. In the B meson rest frame they read 



1 d 3 p x d 3 pi d 3 p u d 3 k 
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64 tt^itib 
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k) \M[ 



Id 



Ei E u Efi 

A^g| 2 + 2^e[A^[]A4^ d (A)] + 2 \M° f " In ^) dE x dE e . (2.14) 



7T 



Integrating eqs. (|2.14|) results in the next-to-leading order total decay rate. Comparing with 
the total tree- level decay rate, the integral over phase-space of eq. (|2.5j) . yields the long-distance 
enhancement factor due to next-to-leading order effects. It is 



r = (1 + ^ + 40 rg = r + r l + r\ + o(a 2 G F ). 



(2.15) 



with 5 s d = — In from eq. (j2. 12)) . The exact form of the next-to- leading order matrix elements 
can be found in App. Q3J 
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2.3. Structure dependent terms 

This section discusses the arising additional electromagnetic next-to-leading order corrections 
that cannot be grasped by simply replacing — > to arrive at a U(l) em gauge invariant 
phenomenological Lagrangian, as was described in the previous section. These include both 
deviations arising from the point-like meson-photon interaction assumed above and additional 
terms arising in the interaction of an off-shell hadronic current. Nonetheless, it is clear that in 
the relevant phase space region for the total inclusive decay rate, namely the region as k — > 
near the infrared divergence, both the real and virtual next to leading order matrix elements 
are completely determined by the leading order decay, and the above procedure accurately 
reproduces the full theory in this region [5,6,25j. Hence, the real emission squared amplitude 
in this limit reads 

\m\\ 2 fc ±° -e 2 f/^-^YK| 2 with PM E{p B , Px }, (2.16) 

depending on whether B or X is charged. Physically this roots in the fact that the wavelength 
of an infinitely soft photon is much larger than the size of any strongly bound hadron. It, thus, 
cannot resolve its substructure and effectively interacts with its summed, then point-like charge. 
Further, such soft photons cannot push the hadronic current significantly off-shell, such that 
off-shell current interactions can become sizeable. 

Introducing non-point-like meson-photon interactions does not only lead to corrections due to 
the hadron's size and its internal charge distribution, it also leads to additional vertices of the 
type X — > X*"f, where X* is a higher resonance of the X meson. This necessarily also introduces 
additional terms in the interaction of the hadronic and the leptonic current, especially if the 
resonance differs in its spin. B* resonances in the initial state occur in the unphysical region. 
Hence, they are only relevant if their width is comparable to, or larger than, the mass separation 
m 2 B — mg, to the initial state B meson. In contrast, D* resonances, for example, occurring in 
a final state line are allowed to be on-shell for a range of photon energies. Thus, a considerable 
correction may arise. [26j find the D* + — > D + j coupling to be compatible with zero while the 
D*° —7- Z)°7 coupling is small, but considerable. Both are considered and discussed in detail in 
Sec. PI 

Generally, following the argumentation of [27, 28] the electromagnetic current of the hadronic 
system can be split into two components: inner-bremsstrahlung (IB) contributions, which ac- 
count for photon radiation from the external charged particles and are completely determined 
by the non-radiative process, and structure- dependent (SD) contributions, which describe inter- 
mediate hadronic states and represent new information with respect to the IB contributions. 
The amplitude of a semileptonic B meson decay with full electromagnetic corrections reads 

A v = i e ^=V xh u u P R ^ ( - ( 7v # + 2p lv ) + V„ u - v l > ( 2 - 17 ) 

with the hadronic current H^, as introduced in eq. (|2.3j) . The hadronic vector and axial form 
factors of the photon-emitting hadronic current, incorporating among others the X —> X*j 
coupling, are given by the unknown non-local operator 

V - V = / d ' x etk ' X ( X \ T MO) Jr(x)} \B) . (2.18) 

J^ m denotes the electromagnetic current, the quark-current in position-space and k is the 
photon momentum. The vector and axial-vector operators of eq. (|2,18p obey the electromagnetic 
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Ward-identities, obtained by contracting k v A v of eq. (|2.17p . 



k v V t 



k u A 



fJlU 



H 
0. 



I-'- - 



(2.19) 



These properties of the individual pieces of the amplitude in conjunction with Low's theorem 
[5 )125] lead to the fact, that the leading terms of next-to-leading order amplitude in powers of the 
photon four- momentum k, i.e. the terms proportional to k~ l and k°, are completely determined 
by the on-shell form factors of the tree-level decay. 

Following [27], corrections beyond O(k ) can be included by separating the non-local opera- 
tor eq. (|2.18p into SD and IB contributions. Since the IB and SD describe different physical 
mechanisms they are separately gauge invariant. Further, the SD amplitude must be of 0{k) or 
higher. This, however, does not prevent the IB amplitude from containing terms of 0{k) and 
higher as well. 

Splitting the amplitude under these restrictions allows more terms to be collected in the IB part, 
still using only the knowledge of the non-radiative matrix element. This offers the advantage to 
obtain more precise predictions for the decay process, without formulating the (mostly unknown) 
SD contributions. The splitting of the transition matrix element requires a corresponding 
splitting of the non-local operator eq. (|2,18p into SD and IB parts. 

The axial contributions are strictly zero for photons emitted from the B- or X-meson, and 
therefore can be considered purely SD. They can be written in the form [2S] 



A 



fj,l> 



A SD 



i e 



fivpcr 



A lP p x k° + A 2 kp( Pi + Pu y 



i tvXpaPx kP (pi + Vvf A z (pi + p v ) + Aip^ 



(2.20) 



Note that the Lorentz-invariant scalars A{ are non-singular in the limit k — > by construction 
and are functions of the three independent scalar variables that can be built with P b, px and 
k. The decomposition of the vector current reads 



V, 



ym + ySD 



where the IB piece is chosen in such a way, that 



k v V, 
k u V, 
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\iv 
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(2.21) 



(2.22) 



Thus, the decay amplitude separates as 



Ay 



G F 

ie-j=V xh u u Pr7 m 



2pi ■ k 



IB 



(2.23) 



Herein, Vj^ can be constructed from leading order information only, cf. App. |Bl The SD vector 
contributions, on the other hand, contain additional information. They can be written as [3D] 



V, 



SD 



kfxPxu ~ (p ■ kjg^ 



+ V 2 



+ v 4 



kfi(pi + Pv)v - (k ■ (pi + Pv))gv,v 

(k ■ (pi+p v ))(p£+Pu)fiPXu - (px ■ k)(pi+ p u )p(pe+ Pu)u 
(k ■ (pi + p v ))(pxpPXv - (px ■ k)p X p(pe + Pu) 



(2.24) 
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where the Lorentz-invariant scalars V{ are functions of the three independent scalar variables 
that can be built with ps, px and k. All IB and SD contributions are finite as k — > 0. The IB and 
SD parametrization of [31J based on [52] can be obtained by a change of basis in eq. (|2.24|) . and 
correspondingly shifting terms of 0(k) and higher into the SD contributions. The IB coupling to 
the electromagnetic current can be used to construct the next-to-leading order matrix elements, 
receiving in principle additional corrections from the SD coupling. 

The knowledge of the full SD contributions for semileptonic B meson decays is modest: [31] 
discusses the matter for B — > irlvj decays, using the soft-collinear effective theory to isolate 
the expressions for the SD contributions in the soft-pion and hard-photon part of phase-space. 
Sec. 14.31 compares their findings with the pure IB result from this study. The differences are 
non sizable. The recent work of [28] addresses the real SD corrections to B — > D i v 7 decays by 
using lattice results of the D* — > D 7 coupling to estimate the dominant SD contributions when 
the D* is on-shell. Sec. 14.31 also compares these SD contribution with the complete SD+IB 



picture. Again the differences turn out to be non sizable. The SD corrections to B — > Dq Iv^f 
are unknown, but given the large widths of the Dq and D\ states a non-negligible correction to 
the pure IB prediction can be expected. 



2.4. Soft-resummation and inclusive exponentiation 

This section discusses a systematic improvement of the fixed order results discussed in the 
previous section. Centring on the exponentiability of soft-radiative corrections and following 
the approach of Yennie, Frautschi and Suura [6], the fully inclusive decay rate 
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i=l \ i=l 

by separating the universal spin-independent infrared divergent terms from the virtual and real 
emission amplitudes. d<& Pf and d<E>fc are the leading order and (multiple) extra emission phase 
space elements, while tir and ny count the additional real and virtual photons present in each 
amplitude. Therefore, using the same convention as before the sub- and superscripts of the 
(squared) matrix elements A4, M and f3 denote their real emission photon multiplicity and 
their order of a in the perturbative expansion relative to the leading order. 
The separation of infrared divergences proceeds by splitting 

Ml 



2^\ M h 2 = s«=)K i2 ■ • 



aBM° Q + Ml , 
S(k)\M°\ 2 + (3i(k), 



(2.27) 
(2.28) 



wherein Mq and f3{(k) are free of any infrared singularities due to virtual or real photon emis- 
sions. This separation can be continued iteratively, leading to 
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(2.29) 
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and 



.2(2vr) 2 
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i=l 
n R 

+ E 



nR 

+E 
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S(h)S(kj) 
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(2.30) 



with /3 nfl = Y^n v =o PnR +nR ■ Note that for a given phase space phase space configuration 
{pi, . . . ,p n , ki, . . . , k nR } the infrared subtracted squared matrix elements /3i(/cj) involve a pro- 
jection onto the single emission subspace {Vpi, ■ ■ ■ , Vp n> Vki}. Of course, momentum conserva- 
tion holds for each projected subset. Thus, for every radiated photon the (3i(ki) are evaluated 
as if this photon was the only one in the event. Hence, truncating the pertubative series in the 
j3 nR at the next-to- leading order leaves every single photon emission correct at 0(a). 
Exponentiating the integral of the eikonal S(k) upon insertion of the identity of eq. (|2,3U|) 
over the unresolved phase space containing the infrared singularity gives rise to the Yennie- 
Frautschi-Suura form factor 



Y(Q) = 2a{B + B{VL)) with 2aB(Q) 



d 3 k ~ 



S(k) 



(2.31) 



n v +n R 
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and the residual perturbative series of the infrared-subtracted squared amplitudes $, 
eq. (|2.26|) . Hence, photon emissions contained in the unresolved soft region O are assumed to 
have a negligible effect on differential distributions, but are included in the overall normalisation. 
Furthermore, Y(Cl) is UV-finite and, thus, does not interfere with renormalisation of the f3™ v+riR 



R 



3. Methods & Implementations 

In this section a short overview over the generators used in this study and their underlying 
principles is given. 



3.1. Blor 

BLOR [18] is a fixed order Monte Carlo Event generator specialised on Qed corrections in 
semileptonic B meson decays. It separately generates decays for the born (1 — > 3) and real 
emission (1 — > 4) phase space events according to their individual cross sections 

dr^ 3 = dr° + dri 

dr^ 4 = dr}. (3.i) 

The infrared divergences present in dFg and dT\ are regulated introducing a small but finite 
photon mass A set to 10 -7 GeV. The expressions of Sec. 12.21 are altered accordingly. 



3.2. Sherpa/Photons++ 

The SHERPA Monte Carlo [19J is a complete event generation frame work for high energy physics 
processes. Although its traditional strengths lie in the perturbative aspects of lepton and hadron 
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colliders it also encompasses several modules for all non-perturbative aspects. In this work the 
hadron decay module HADRONS++ [33] an d the universal higher order Qed correction tool 
PHOTONS++ [2D] are used. 

HADRONS++ generates the leading order decay events according to the respective form factor 
parametrisations of the involved hadronic and leptonic currents. PHOTONS++ then rewrites the 
differential all orders inclusive decay width of eq. (|2.26p as 



r - r. E ^ / »■ *(*) ft [SIW*. «)] ( i + 1 + g + ^ 



A) i=l A) ' 



(3.2) 



cf. [20]. Hence, the leading order input is corrected both for soft photon effects to all orders 
and hard photon emission to an arbitrary order. Jp < 1 represents the various Jacobians 
occurring when factoring out the leading order term. The perturbative series in the infrared- 
subtracted squared matrix elements includes only terms up to 0(a) for this study. The Nlo real 
emission squared amplitudes can alternatively be approximated using Catani-Seymour splitting 
function 



i<j 

where in i and j run over all particles in the process. The exact form of the can be found 
in [20]. This approximation is also used if the exact real emission matrix element is not known. 
The infrared cut-off was set to 10 _6 GeV in the rest frame of the charged dipole, i.e. the B + — £ + 
or the X~— i + system, respectively. 

3.3. Photos 

The PHOTOS Monte Carlo [15] is an "after-burner" algorithm, which adds approximate brems- 
strahlung corrections to leading order events produced by an external code. In this study it is 
taken as a reference for the differential distributions since it is widely in use. PHOTOS bases on 
the factorisation of the real emission matrix element in the collinear limit 

l^i I' = EMI'" fM. (3.4) 

i 

The radiation function f(pi,k) is given to leading logarithmic accuracy. It incorporates the 
Altarelli-Parisi emission kernel for radiation off the particular final state particle. Its exact form 
is spin dependent and can be found in [15pi6j . In its exponentiated mode the number of photons 
follows a Poissonian distribution, while the individual photon's kinematics are determined by 
applying above equation iteratively. To also recover the soft limit of real photon emission matrix 
elements an additional weight was introduced [37] 



W soit 



nc * 



n~{ nc 

i=ij=i 



Pi- 



p r ki 



(3.5) 



wherein n 7 and uq are the photon and the final state charged multiplicity of the process. To 
regularise the emission function in the soft limit an energy cut-off is imposed in the rest frame 
of the decaying particle. The collinear divergence is regularised by the emitter's mass. 
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Initial state radiation is not accounted for as the mass of the decaying particle is the largest scale 
in the process and there are no associated collinear divergences or logarithmic enhancements. 
It could only be accounted for by supplementing PHOTOS with a matrix element correction. In 
case of heavy initial states eq. (j3.5[) still approximately recovers the soft limit. 
In this analysis PHOTOS version 2.13 has been used in its exponentiated mode, including the soft 
interference terms. The infrared cut-off was set to W~ 7 m^o,+ , respectively. BLOR supplemented 
the leading order decay events. 



4. Results 



Parameter Value 



m T(4S ) 

r T(45) 

m B o 
m B + 
m D o 
m D - 
m n *o 
m n -o 
m n o 



10.5794 GeV 
20.5 MeV 
5.27950 GeV 
5.27913 GeV 
1.86950 GeV 
1.86484 GeV 
2.403 GeV 
2.352 GeV 
0.134976 GeV 
0.13957 GeV 



Parameter Value 



m e 

mw 
a 

G F 



0.0005109989 GeV 
0.10565837 GeV 



80.419 GeV 
1/137.035999 
1.16637 • 10~ 5 GeV" 



Tab. 



1 Parameters used for all inclusive and differential decay rate calculations. All particle widths, 
except the T(45) width, are considered negligible. 



In this section the results of the inclusion of the complete next-to- leading order corrections, 
both real and virtual short and long distance contributions, are reviewed for different decay 
channels with different form factor parametrisations of the hadronic current. First, results 
for the next-to- leading order inclusive decay rates and their effects on the extraction of |Vxb| 
from decay measurements will be shown in Sec. 14.11 Then, in Sec. 14.21 the effects on differential 
distributions and spectra are investigated and compared against the standard tool used in many 
experimental analyses, PHOTOS. The parameters used are detailed in Tab. [TJ 



4.1. Next-to-leading order corrections to decay rates 

One key prediction of this publication is the process specific correction factor for higher order 
electroweak effects. These higher order corrections enter measurements of the Ckm mixing 
angles via 

Pmeasured — Vqcd Vew I ^xb | PlO j 

and, thus, 



|Kb| = — • J " • ( 4 -2) 

^ew y r^ CD r LO 

Plo is the leading order phenomenological decay rate stripped of the Ckm mixing angle. ?7qcd 
and r; EW incorporate the higher order Qcd and electroweak corrections. Both contributions 
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B°- 

B+ 

B+ 



B°- 

B+ 

B+ 



B°- 

B+ 

B+ 



D~ e+ v e (7) 

D- V + u„ (7) 
► D° e+ v e (7) 



1.0222 
1.0222 
1.0146 
1.0147 



7r e+ v e (7) 

7T~ fi + Vjx (7) 

► 7T° e+ ^ e (7) 

> 7T° /i+ ^ (7) 



1.0555 
1.0545 
1.0411 
1.0401 



' Vfi (7) 

£>SV + ^(7) 



D* 
D* 



*o 
' 
s*o 



1.0224 
1.0226 
1.0142 
1.0144 



l/%w 

1±2±17± 1) 0.9891(1 ±1±4±1) 
1±2±17±1) 0.9891(1 ± 1 ±4± 1) 
1±1±39±16) 0.9928(1 ± 1 ± 10 ±4) 
1±1±39±16) 0.9927(1 ± 1 ± 10 ±4) 



1±4±148±48) 0.9734(1 ± 1 ± 33 ± 10) 

1±4±136±48) 0.9738(1 ± 1 ±31 ± 11) 

1±3±100) 0.9801(1 ± 1 ±23) 

1±3±89) 0.9805(1 ± 1 ±21) 



1 ± 2± 
1 ±2± 
1 ± 1 ± 
1 ± 1 ± 



10) 0.9890(1 ± 1 ± 2) 

10) 0.9889(1 ± 1 ± 2) 

35) 0.9930(1 ± 1 ± 8) 

35) 0.9929(1 ± 1 ± 8) 



Tab. 2 Predictions for rj EW = y/1 + i5 s d + 5id for the summed next-to-leading corrections are listed. 

The uncertainties in the parentheses are the sum of numerical, next-to-next-to-leading order, 
matching and missing structure dependent contributions. For B + — > ir° I v (7) and B°' + — > 
Dq '° I v (7) decays no predictions for structure dependent contributions are known. 



factorise at the Nlo level. The electroweak correction factor is determined in this study at 
Nlo accuracy in the QED-improved phenomenological long-distance description of the hadronic 
decay and at leading logarithmic accuracy in the underlying short-distance partonic decay in 
the Standard Model. While the leading logarithm of the short-distance correction depends only 
on the matching scale A of both descriptions and is, thus, independent of the actual decay 
properties, the long-distance corrections, in contrast, are sensitive to exactly these specifics. 
The resulting correction factors, 

rL = l±4d±5i d = l + ^^ + 0(a 2 ) (4.3) 

1 

are a central outcome of this study. They are presented in Tab. [2] for the different semileptonic 
decay channels of charged and neutral B mesons into (pseudo) scalar mesons. As is evident, 
the long-distance Qed corrections break the strong isospin symmetry. This originates in the 
different masses of the charged mesons in the strong isospin rotated decays. These masses both 
determine the amount of radiation from the meson line and enter the loop integrals. For the 
same reasons the corrections for pions are larger than for D mesons. On the other hand, the 
difference in the size of the correction between the two leptonic channels of each decay mode is 
only very small since both masses are insignificant compared to the hadronic mass scales. 
The total uncertainty Ototal of the summed long- and short-distance correction is given by 

a total = ^numerical + °"nnlo ± °A ± ^SD ■ (4-4) 

The leading uncertainty originates from the matching of the short- and long-distance results. 
Here, the mismatch of both theories, the Standard Model for the short-distance corrections and 
the QED-enhanced effective theory for the long-distance corrections, at the scale A, as discussed 
in Sec. I2.2.1| is the main source. These matching uncertainties are estimated choosing A of 
the scale of the final state meson's mass as a central value for A and then taking the difference 
to the result when using the scale 2A. Note that half the charged final state meson's mass 
is not a sensible choice for the matching scale since in this case a considerable fraction of the 
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n 2 
VEW 


1 / 

l/f?EW 


£> u — i 


► D Iv 


1.0222(17) 


0.9891(4) 


B+ - 


> D°£u 


1.0146(43) 


0.9928(10) 


B 




1.0186(29) 


0.9909(7) 






^EW 


V^ew 




7T~ £ V 


1.0550(150) 


0.9736(34) 



Tab. 3 Averaged integration results for r;| w = 1 + <5 s d + #id and l/r?Ew : The uncertainties in the 
parentheses are the sum of numerical, next-to-next-to-leading order, matching and missing 
structure dependent contributions. 



Measurement 


£(1) \Vch\ x 10~ 3 with 
sd corrections only 


£(1) l^bl x 10~ 3 with 
sd and Id corrections 


BABAR tagged [38] 


42.30(2.36) 


42.21(2.34) 


Measurement 


Vub x 10~ 3 with 
no corrections 


F u b x 10~ 3 with 
sd and Id corrections 


BABAR untagged [39] 


3.60(±°;g) 


3.50(1°;^) 



Tab. 4 The impact of long-distance corrections on the correction factor ?y EW is shown. Both measure- 
ments estimate real next-to-leading order corrections via PHOTOS and, thus, use phase-space 
cuts and variables that provide a reduced sensitivity on the modeling of final state radiation. 
Q{\) is the normalisation of the heavy quark form factor. Note that |Kib| aimed analyses 
usually do not apply short-distance corrections and that the form factor prediction of [40] 
was used to extract \V U \\ from the measured B a — > 7r~ £ + v% partial branching fraction. 

real radiation cross section would occur at scales greater than A and being described by the 
phenomenological model, thus ending up in the wrong picture. 

Along these lines, also again referring to Sec. 12.2.11 the results for decays into pions should 
be considered with care: apart from the afore mentioned conceptual problems introduced by 
the large separation of both hadronic mass scales, large parts of the real emission phase space, 
described in the phenomenological model with a point-like pion, are actually able to resolve 
the pion, the charged and the neutral one. Thus, the model used in this study is not entirely 
accurate for such decays. 

Further, the effect of additional real short-distance contributions is studied by comparing the 
pure inner-bremsstrahlung calculation with the real emission results of [31 J and [28] • This, 
however, can only be done for and B — > D Iv (7) and B° — > tt~~ t + v (7) decays. Such terms are 
unknown for the other decay modes. Consequently, no error can be associated with them. Fi- 
nally, the electroweak next-to-next-to-leading order effects are estimated as cr nn i = a (<5 s d + 5id)- 
It is clear that the complete model dependence from our phenomenological treatment of the 
long-distance contributions cannot be grasped with the above error estimation. 
Tab. [3] presents the same results for B — > D£v(^) and B° — > tt~ i + v (7) averaged over the 
different lepton species and the isospin rotated decays. The isospin averaged result is corrected 
for the difference in the production rate of B° and B + mesons, i.e. 

<5 sd + fid = Osd + S+) /+- + s d + 5i° d ) foo , (4.5) 
with &t and 5® d the isospin breaking contributions for charged and uncharged charmed decays 
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extracted from Tab. [2J The latter correction factor is the adequate choice to correct |V^b| from 
measurements demanding isospin. 

The impact of these correction factors on two selected measurements are listed in Tab. [H Note 
that both measurements used PHOTOS to estimate the effect on radiative corrections. Applying 
the stated factors only corrects the overall normalisation, differences due to changes in kinematic 
distributions, see Sec. 14.21 result in another correction for the extracted values of |Vdb| an d |Kib| 
that cannot be estimated here. The correction for \Vub\ shown here is therefore an illustration 
only and should be used with care. 

4.2. Next-to-leading order corrections to differential rates 

In this section the results of both BLOR and SHERPA/PHOTONS++ are presented and compared 
against PHOTOS. The focus lies on the absolute value of the spatial momentum of the produced 
meson and the lepton, i.e. \px\ and \pe\- The chosen frame for these observables is the centre-of- 
momentum system of the electron and positron beam. Thus, the decaying B° and B + mesons 
already carry momentum corresponding to the Y(4S) — > B°B° and Y(4S) — > B + B~ decay 
kinematics. All quantities are shown as bare quantities, i.e. no recombination of photon and 
lepton/meson momenta was used. This is directly applicable if the charged particle momenta 
are extracted by measuring their curved tracks in a magnetic field or in photon- free calorimeters, 
as is the case in many BABAR and BELLE analyses. 

The prediction of each generator is normed to its inclusive decay width and the ratio plots show 
the relative difference 

i dr\ _ l dr rcf 

a r to t,i dO rtot.ref dO 

° ~ i dr, i dr rcf 

Ttot,i dO rtot.ref dO 

to the PHOTOS prediction for the given observable O in the given bin. Hence, the short- 
comings of the approximations inherent in the standard tool currently used by most experiments 
are plainly visible. Further, due to the choice of normalisation, systematic errors, shown to 
be dominant in the previous section, are negligible here. Hence, the error bands shown are 
statistical errors only and are of comparable magnitude for all three generators predictions. 

4.2.1. Decays B -> Biv t 

In Fig. [5] the predictions of all three generators for the decay B° — > D~ e + v e using a form factor 
parametrisation from Heavy Quark Effective Theory, cf. App. IA.ll are presented. For most of 
the phase space the agreement of the next-to-leading order shape is good. However, as can be 
seen in the absolute lepton momentum plot, SHERPA/PHOTONS++ predicts a slightly different 
shape close to = 0. This limit cannot be measured directly in most experiments. 

Nonetheless, it influences the extrapolation to the full phase space, and thus determinations 
of total decay widths. In [JT], for example, uncertainties related to such an extrapolation 
constituted up to 50% of the total experimental error. Further, the slight slope of the order 
of 1% exhibited in the prediction of SHERPA/PHOTONS++ for \px\ is also present in the \pg\ 
spectrum. Both differences originate in the different modeling of radiative energy loss between 
the different programs, as shown in the following. 

Fig. [6] displays the radiative energy loss, i.e. the sum of the energies of all photons radiated, in the 
rest frame of the decaying B°. The predictions of the three different generators are shown in the 
left panel. For single photon emission in such a final state dipole there exists a kinematic limit 
on the photon's energy: E™ x = rn x +m e \ ^ the dipole > s rest frame or gmax = 
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I in the e + e rest frame 




\p(\ in the e + e rest frame 



Fig. 5 




Lepton and Meson momentum spectrum in the e + e~ rest frame in the decay B — > D~ e + v e . 
All spectra are normed to the total inclusive decay width predicted by the respective generator. 
The ratio plot gives the relative deviation, bin by bin, of the predicted shapes with PHOTOS as 
reference. The shaded yellow area gives the statistical uncertainty of the reference distribution. 



£ 7 in the B rest frame Ey in the B rest frame 




E 7 [GeV] E 7 [GeV] 



Fig. 6 Total radiative energy loss, i.e. the sum of all photons radiated, in the decay B° —> D~ e + v e 
in the B rest frame. All spectra are normed to the total inclusive decay width predicted by 
the respective generator. The left panel shows the predictions of all three generators and 
the PHOTOS prediction is taken as the reference in the ratio plot. The right panel shows the 
predictions of SHERPA/PHOTONS++ in its full YFS<8>Nlo exponentiated mode (green), a mode 
where the exact Nlo matrix element of the perturbative expansion is replaced by universal 
Catani-Seymour dipole splitting kernels (red) and two modes where the exact real emission 
matrix elements are used, but the expansion in the resolved emission region is truncated at 
0(a) (blue) and 0(a 2 ) (orange), thus allowing at most one and two photons, respectively. 
Here, the full exponentiated SHERPA/PHOTONS++ prediction is taken as reference in the ratio 
plot. 
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2.3083GeV in the B rest frame. This limit is clearly visible in Fig. All events exceeding it, 
i.e. radiating more than -E™ ax , must exhibit multi-photon radiation. Here (at least) two hard 
photons recoil against each other. Hence, this feature is present in both the SHERPA/PHOTONS++ 
and PHOTOS predictions, but not in the fixed order Nlo prediction of BLOR. This tail is an 
0(a 2 ) effect in the hard radiation and can only be described approximatively here. In PHOTOS 
it is described by an iteration of the emission kernels while SHERPA/PHOTONS++ describes this 
part of the spectrum by the next-to-leading order hard emission amplitude j3\ summed over all 
projections onto the single emission subspaces, cf. Sec. 12.41 and 13.21 

On the other hand, multi-photon radiation also enhances the amount of radiation in the region 
Ely < £™ ax , if the probability of two relatively hard photons is sufficiently large. As exemplified 
in Fig. in B° — > D~ e + u e double photon emission is relatively probable and, hence, leads to 
such an enhancement, whereas due to the much larger muon mass this feature is nearly absent 
in the decay to muons, Fig. [9] Of course, this enhancement of the radiative energy loss is an 
effect of 0(a 2 ) and can therefore only be described approximatively here. 
Of equal importance as multi-photon radiation is the presence of the exact real emission matrix 
element, as is also shown in the right panel of Figs. [6] and [9] Approximating the real emission 
matrix elements with Catani-Seymour splitting functions, reproducing the Altarelli-Parisi split- 
ting functions used in PHOTOS in the (quasi-) collinear limit, leads to a mis-estimation of the 
radiative energy loss in the regime close to the kinematic boundary. It seems, however, that in 
the present cases collinearly approximated multi-photon emission mimics the exact fixed-order 
Nlo behaviour reasonably well. 

Close to the kinematic boundary on single photon emission, E max , the vertex emission diagrams 
become important, as do the corrections for t ^ t' (cf. App. [B]). These corrections have different 
sizes for the electron and muon channels due to their different masses and radiative properties. 
In principle, here also the structure-dependent corrections of Sec. 12.31 play a role. But, as is 
investigated in Sec. l4.3l thev have negligible impact on the shape of the differential distributions. 
Thus, they can be safely neglected here. 

Further, both SHERPA/PHOTONS++ and PHOTOS share a common soft limit, showing the com- 
patibility of the inherent soft resummation of SHERPA/PHOTONS++ and the superimposed soft 
limit correction in PHOTOS. BLOR exhibits an (almost) constant off-set of a few percent owing 
to the lack of resummed contributions modeled in SHERPA/PHOTONS++ by virtue of the YFS 
form factor. 

The lepton and meson spectra for B° — > D~ /i + are shown in Fig. [8] Figs. [10] and [11] show 
strong isospin rotated processes B + — > D° e + v e and B + — > D° fi + f M , respectively. Thus, the 
radiating dipole is spanned between the initial state B + and the lepton. Radiation off the initial 
state meson is suppressed by its much larger mass, as compared to the D~ . Thus, multi-photon 
emission is also strongly suppressed. In these cases SHERPA/PHOTONS++ predict slightly smaller 
radiative corrections in the electron decay channel than either PHOTOS or BLOR. The differences 
are of the order of five percent; note that the scale was enlarged in the reference plot to better 
highlith the differences. These, to the largest extent, root in differences in the modeling of 
emission off the initial state charged meson. 

4.2.2. Decays B ->• D%iv t 

The final state lepton and meson momentum spectra in the decays B°' + — > Dq~'° £ + U£, with 
£ = e,fx, are shown in Figs I12L I13|. and[T4"l [T5l The B —> Dq transition current is modeled using 
Leibovich-Ligeti-Stewart- Wise-parametrised form factors, cf. App. IA.31 Except for differences 
in the form factor parametrisations due to the Dq meson being a scalar instead of a pseudo- 
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photon multiplicity with Ej > 30MeV in the e + e rest frame 



:° tt 



Photos 

Blor 

Sherpa 



photon multiplicity with E } , > 30MeV in the e + e rest frame 



o.X 
0.6 
0.4 

0.2 
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Blor 

Sherpa 



Fig. 7 Multiplicities of photons with at least 30MeV in in the e + e~ rest frame B° — > D~ e + v e on 
the right hand side and B° — » D~ /i + on the left hand side. In the ratio plot PHOTOS was 
chosen as the reference. 
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Fig. 8 Lepton and Meson momentum spectrum in the e + e rest frame in the decay B° — > D 
The PHOTOS prediction is taken as the reference in the ratio plot. 
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£ 7 in the B rest frame 



£ 7 in the B rest frame 
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E y [GeV] 

Fig. 9 Total radiative energy loss, i.e. the sum of all photons radiated, in the decay B° — » D 
in the B rest frame. The labels are identical to those in Fig. [51 
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\px \ in the e + e rest frame 
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Lepton and Meson momentum spectrum in the e + e~ rest frame in the decay B —, 
Dq ~ n + v^. The PHOTOS prediction is taken as the reference in the ratio plot. 




24 



\px\ in the e + e rest frame 



£P 0.14 




in the e + e rest frame 



0.04 



o 

0.04 
0.03 

0.02 
0.0 1 



-0.01 

-0.02 

-0.03 §=- 
-0.04 



o 



LO 

Photos I 
Blor I 
Sherpat 



0.5 



i-5 



2-5 
|p,| [GeV] 



Fig. 15 Lepton and Meson momentum spectrum in the e + e rest frame in the decay B + 
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The PHOTOS prediction is taken as the reference in the ratio plot. 



scalar and its higher mass, the effects of higher order corrections are comparable to the case of 
Sec. ELD 

4.2.3. Decays B -> ir£ve 

Figs. [To] [T71 and {TE[ [T9l show the decay channels into charged and neutral pions. The B — > n 
transition current is modeled using the Ball-Zwicky form factor model described in App. IA.21 
Here, because of the comparably small mass of the charged pion effects due IB corrections for 
t 7^ t' become important. The structure-dependent contributions still have negligible impact on 
the differential distributions, as is shown in Sec. 14.31 Again, in the electron channel of the decay 
into a charged pion the correct treatment of hard multi-photon radiation, assuming they are 
sufficiently well described in the QED-enhanced phenomenological model, leads to comparably 
large deviations. 

Nonetheless, it should be noted, as was also discussed earlier, that the matching procedure 
employed in this study runs into conceptual problems when applied to a B — > tt transition due 
to large difference between the hadronic mass scale A = and the maximal energy of an 
emitted photon, £™ ax = 2.6379GeV (B° — > TT~e + v e ). Consequently, a considerable fraction of 
the real emission phase space wherein the photon is able to resolve the pion is described by 
the effective theory only. Thus, the results obtained here should be considered with caution. 
However, they still are an improvement over the leading logarithmic corrections employed in 
standard analyses. 

4.3. Influence of explicit short-distance terms 

In section the influence of explicitly calculated structure dependent terms, as introduced in 
Sec. 12. 3[ is investigated. The analysis is performed for the decay channel B° — > ir~ e + v e with 
the results of [31] and the decay channels B 0,+ — > D~'° e + v e with the results of [28]. The same 
conclusions also apply to the muon channels. 

For the charged pion channel, [31] employs a form factor model of the heavy-hadron chiral 
perturbation theory (HHxpt) [421-145] . valid in the region of E~ > lGeV. Despite the mismatch 
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|px| in the e + e rest frame \p{\ in the e + e rest frame 




o 0.5 1 1.5 2 2.5 3 o 0.5 1 1.5 2 2.5 3 

Ipxl [GeV] \f t \ [GeV] 

Fig. 16 Lepton and Meson momentum spectrum in the e + e~ rest frame in the decay B° —> ir~ e + v e . 

Both matrix-element-corrected multi-photon radiation and the IB terms for t 7^ t' exhibit a 
strong influence here. The PHOTOS prediction is taken as the reference in the ratio plot. 
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Fig. 17 Lepton and Meson momentum spectrum in the e + e~ rest frame in the decay B° — > ir~ /i + v^. 

The IB terms for t 7^ t' exhibit a strong influence here. The PHOTOS prediction is taken as 
the reference in the ratio plot. 



£sd + £ld(IB) 5 s d + ^ld(IB+part.SD) °"SD 

5 U 1/^(7) 0.02223(6) 0.02225(7) 0.00002 
ff+ L>° 1+ ^(7) 0.01463(5) 0.01627(6) 0.00158 

Tab. 5 The effect of including the partial SD terms arising due to intermediate excited D*~ and 
D*° mesons is shown. Because of the unknown size of the full SD contributions this single 
term, argued to be dominating, is used as an estimate on the systematic uncertainty ergD 
associated to the IB-only result of Sec. 14.11 
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\px\ in the e + e rest frame \pg\ in the e + e rest frame 
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Fig. 19 Lepton and Meson momentum spectrum in the e + e~ rest frame in the decay B + — > tt° /i + v^. 
The PHOTOS prediction is taken as the reference in the ratio plot. 
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Fig. 20 The photon energy spectrum in the decay B° — > ir~ e + v e 7 (E 7 > lGeV) is shown. The 
complete IB+SD result of [21] (blue) is compared against the prediction of the IB terms only 
according to Sec. 12.31 (black) in the HHxpt form factor model. In the ratio plot the latter is 
taken as the reference. 



£ 7 in the B rest frame Ey in the B rest frame 
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Fig. 21 The photon energy spectra in the decays B° — >• Z)~ e + i/ e 7 on the left and B + — ► i? e + v e 7 
on the right are shown. The result including partial SD terms arising due to intermediate 
excited mesons, D*~ and D*° [55], and the associated coupling X* — > Xj (grey) is compared 
against the prediction of the IB terms only according to Sec. 12.31 (black) in the HQET form 
factor model. In the ratio plot the latter is taken as the reference. 
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of the form factor model used in the present study the result depicted in Fig. [201 shows the 
structure-dependent contributions have little influence on the photon energy spectrum in this 
region. Due to their finiteness in the limit k — > 0, they are expected to behave similarly for 
Ely < lGeV. For the Ball-Zwicky form factor model used in this study, cf. App. IA.21 the IB 
and SD correction are expected to behave similarly. 

In the D meson channels, [28] uses lattice results for the trilinear couplings of an exited D* 
meson to a photon and its ground state: 5d*+d+ 7 = —0.1(7) and to* _D 7 = 2.7(1.2). The 
effects manifest themselves as corrections to the total decay widths and are summarized in 
Tab. [5j they prove minor in the case of D~ , and sizable in the case of D°. The case is similar 
for the radiative spectra: a slight change in the shape of the radiative energy loss for the D° 
channel on the scale of less than 5%, while no such change occurs in the D~ channel. These 
particular SD corrections, however, only form a single term in one class of SD correction. 
Note that higher order charm resonances, i.e. through D**D~f processes, do not contribute to 
real corrections at Nlo of the studied decay modes because of angular momentum and spin 
conservation. Therefore, the lowest order are believed to be one of the dominant terms, they 
are taken as an estimate for the error associated to neglecting all SD contributions. 
Note, that the samples containing (partial) SD terms have larger statistical uncertainties. 

5. Conclusions 

In this paper, electroweak correction to semileptonic B decays were studied. A long-distance cal- 
culation in the Qed improved effective Lagrangian was matched to the partonic short-distance 
result of [8] for (pseudo) scalar final state mesons. Structure dependent terms, e.g. due to non- 
local photon-charged meson interactions, intermediate resonant meson propagators or modifi- 
cations to the effective weak meson decay due to off-shell currents, were not taken into account 
for the computed central values. This was done because they are only known for a very limited 
set of processes, and there usually only partially. 

The results achieved with this method, detailed in Sec. [U give more reliable predictions for total 
and differential decay rates, accompanied by quantifiable errors. The improved predictions of 
the total decay rate were applied on two selected measurements of [38J and [39] resulting 
in small corrections to their respective central values. These exemplifications, however, are a 
mere reweighting of their stated results and correction factors. To fully assess the impact 
of the corrections to the leading order decay presented in this paper the form factors of the 
phenomenological models will have to be refitted with the results presented in Sec. 14.21 Here, 
special attention is again drawn to the large deviations near the endpoints of the kinematic 
distributions arising when both fixed next-to-leading order and resummed leading-logarithmic 
calculations are combined. Finally, the parts of the analyses relying on Monte Carlo estimates 
of the radiation pattern need to be corrected for the improved description presented in this 
publication. 

It should again be emphasised that the results presented for the decay of a B meson into a 
pion should be considered with care. The prescription of matching long- and short-distance 
corrections runs into conceptual problems for this particular process. This is due to the large 
hierarchy of the scales of the pion mass (the scale where a photon is able to resolve a pion) 
and the maximally allowed photon energy (£'™ ax ~ 2.5GeV, cf. Sec. I4.2.3|) . Consequently, the 
prediction for the total decay rate may receive significant corrections when a more elaborate 
matching scheme is used. Nonetheless, the differential distributions are unaffected. 
Similarly, structure-dependent contributions, where known, have been shown to have negligible 
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influence on the differential distributions while their impact on the total decay rates can be 
sizeable. Despite this fact, these structure-dependent contributions were not included in the 
predictions of the central values of the total decay rates, but only have been used to estimate 
their potential error. This treatment is justified since in all cases considered here they are only 
known partially. 
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All histograms were plotted using tools from Rivet [46] . 

A. Form factor models of exclusive semileptonic B meson decays 
reviewed. 

A.l. Form factors for B — > Dtv 



Tab. 6 Parameter values used for the transition current (-D|T^|i?) taken from Heavy Quark Effective 



Parameter Value 



G{1) 0.98 
p 2 D 1.19 



Theory. 



The vector current describing the semileptonic B — >• D tv decay is given by 



(D\V\B) = ^m B m D (h+(w)(v B + v D )^ + h.(w){v B - vof) 
with the heavy quark form factors h± parametrised [17] 



(A.l) 



h+ (w) = 0(1) x [1 - 8p 2 D z + (51p 2 D - 10) z 2 - (252p 2 D - 84) z 3 
h_{w) = 0. 



)] 



(A.2) 
(A.3) 



at w = 1. The values used are given in Tab. [H 




>£, the form factor slope, and G(l) the normalisation 



A.2. Form factors for B — > n £ v 



The vector current describing the semileptonic B — > ir £ v decay is given by 



{pb-p-kYJ fo(t) , 
(A.4) 
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Par ameter Value Parameter Value 

~2 oe /in 7Tn\?2 



mj +i 28.40 GeV^ rf+1 0.744 

40.73 GeV 2 -0.486 



2 



U2 wj./owv r /+a 



m 

nijj" 33.81 GeV 2 r /o 0.258 



Tab. 7 Parameter values used for the transition current (wlV^lB) taken from the pole parametrisa- 
tion in 1301. 



with form factors parametrised as 

r /+l , r f+2 



/+(*) 

/o(t) 



1 



/+1 



1 



— *5 

/+2 



1 



(A.5) 
(A.6) 



fa 



r / + D r f +2 > an d f/o are normalisations and rrif +2 , and m/ pole masses. Their values are 

listed in Tab. 

A. 3. Form factors \ox B D* n tv 



Parameter Value 

0.3571 GeV" 1 
0.1042 GeV" 1 



Parameter Value 



(b 

£ 



-1.0 



A 

A* 

C(i) 



0.4 GeV 
0.75 GeV 
1.0 



Tab. 8 Parameter values used for the transition current (DqIA^IB) taken from the pole parametri- 
sation in 13511351. 



The axial- vector current describing the semileptonic B — > Dq £ v decay is given by 

<£> *K|£) = ^/m B m Dl) (g+ (w) (v B + v D * Q ) M + <?_ (w) (v B - v D * ) M ) . 
with the form factors g± parametrised as [38,49] 

wA* - A] 



5+H 



with 

CiH 



2(^-i)dH-3CH 



w + l 



- e b 



A*(2w + 1) - A(w + 2) 
w + 1 



CH-2(u;-l)CiH 



C(l)x [l + C'( 
A CM- 



1)] . 



(A.7) 



(A.8) 
(A.9) 

(A.10) 
(A.11) 



where £' denotes the form factor slope. The parameters are defined as e c = e& = 2m~> 
A = mo — m c , and A* = mrj* — m c . Their values are listed in Tab. [8) As can be seen from the 
hadronic current of eq. (]A.7p the role of vector and axial- vector terms are reversed in decays to 
scalars as opposed to decays to pseudo-scalars. Thus, in the discussion of IB and SD terms the 
role of V^v and A^ v are reversed. In particular, = and A 1 ® ^ 0. 
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B. Next-to-leading order matrix elements 



This appendix presents details on the real emission matrix elements with special focus on the 
inner bremsstrahlungs (IB) vertex emission terms. The corresponding virtual emission matrix 
elements can be found in |18j . Sec. 4.3. They are calculated in D = 4 — 2e dimensions and 
are UV regularised using the Pauli-Villars prescription [22] of introducing an unphysical heavy 
photon of mass A, the matching scale to the short distance result. 

The summed matrix element of the Feynman graphs a to c in Fig. [3] for B + — > X° i + ve 7 is 
V 2 Pi m k 



i e — I Kb PB € , H^{p B - k, px;t') u u Pr7 m vt 
V2 Pb -k 



1 



I KA 2 (y> 1 \ 

' 1, vertex emission ' v ' / 

Similarly, the summed matrix element of the Feynman graphs a to c in Fig. U] for B° — > 

X~ l + vi 7 is 

M i = ie —m V y^H lJL {pB,Px]t) u v P K ^ ^ V£ 

v 2 Pe ■ fe 



ie^=V xh PX 6 H^(p B ,p x + k;t') u u P R ^ ve 
V2 Px ' k 



1 



"^^1, vertex emission ' C^"^) 

The real emission changes the definition of the four momentum transfer squared, depending on 
the emission leg: 

t = {pb-Px) 2 = (pe+Pu + k) 2 , t' = (pB-px-k) 2 = (pe+Pu) 2 , (B.3) 

and it is 

HM'P2-,t) = {px+P2) lx U(t) + [p l - P2 ) ix /_(t), (B.4) 

as defined in eq. (|2.3|) . The emission terms off the external mesons and leptons in eqs. (jB.ip 
and ()B.2j) are, however, not gauge invariant by themselves. The vertex emission terms are thus 
needed to restore gauge invariance. Assuming t = t' , f±(t) = f±(t'), the vertex emission terms 
of the constant-form-factor Qed invariant Lagrangian of eq. (|2.13p are recovered: 

-^tertex emission = -» e (/+(*) + /_(*)) U U P R ^V£ (B+ ^ X° (B.5) 

V2' 

This is a reasonable approximation for heavy meson processes. Here, hard photons are emitted 
predominantly collinear to the charged lepton. Light meson processes, however, also radiate a 
considerable fraction of their hard radiation in the direction of the light meson and, hence, give 
rise to non-negligible corrections for t 7^ t'. There are two ways to obtain these corrections. 
A: Supposing the form factors f± can be expanded around t, the hadronic current reads 



MU emission = + ™ "4^b (/+(*) " /-(*)) U v P R {*V t (B° -+ X~ £ + U e ) . (B.6) 



fc=o 2 dt' 2 



k=0 
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with t' = t + k' and k! = — 2k- (ps — Px)- Introducing eq. (|B.7|) into eqs. (jB.ip and (jB.2[) and 
employing Ward's identity [50] to obtain the gauge restoring terms, results in 

1 B+^x°e+u e 



Ml 



vertex emission 

i e 



^§Kb u v P R7 % K - e*) (b\ (f + (t) + /_(*)) - 2( PB - VxT^ fc=Q 

0(fc 2 ) , (B 



and 



1, vertex emission 



fc=0 



+ 0(& 2 ), (B.9) 

respectively. Neglecting higher order terms in the expansion of eq. (jB.Tj) results in Low's matrix 
element [5] for these processes. This approach, by implying the existence of a Taylor-series 
representation of the form factors f±(t), yields a consistent result both for the interaction terms 
of the phenomenological Lagrangian and the Feynman rules. However, the exact functional 
form of the form factors has to be known. Further, by including higher order corrections in 
k structure dependent contributions to the matrix element are introduced: the isolated terms 
that restore gauge invariance are not unique, and undesired ambiguities are apparent. The 
impact of such terms were studied for decays in [32,51j, finding negligible impact on the 
next-to-leading order decay rate. This result however can't be extrapolated to B meson decays, 
due to the wide range of possible excited intermediate states. 

B: A result independent of the functional form of the form factors, and thus not relying on their 
differentiability, can be derived similarly. Instead of their argument, the form factors themselves 
are decomposed 

/±(0 = f±(t) + Z±(t,t') . (B.10) 

Now, the missing terms for achieving gauge invariance of the real emission amplitude are de- 
termined as 

1, vertex emission 

= ie^V xh UuPk-fvt (-(/+(*') + /-(*')) e 



V2 



+(PB +Px)» £*• Z + (t,t') + (pb-Px)h e* MM')) (B.H) 



and 



1, vertex emission 



= *e%ib UvPRrfvt ((/+(0 -/-(O) e* 

+(PB+Px)„ €*■ Z+(t,t') + (pb-Px)h e*- , (B.12) 

respectively. Z±(i,f) = ^Z ± (f,f) = k a ^{f±(t>) - f±(t)) = k a Z£(t,t'), k ■ n £ but 
otherwise arbitrary. Through the definition of n, however, again ambiguities are introduced 
into this generic result which, again, are assumed to be negligible in this paper. The emission 
terms constructed this way constitute the inner-bremstrahlungs part of non-local emission term 
V™ of Sec. 
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